Programming scale-free optics in disordered ferroelectrics 
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Using the history-dependence of a dipolar glass hosted in a compositionally-disordered lithium-enriched 
potassium-tantalate-niobate (KTN:Li) crystal, we demonstrate scale-free optical propagation at tunable 
temperatures. The operating equilibration temperature is determined by previous crystal spiralling in the 
temperature/cooling-rate phase-space. © 2012 Optical Society of America 
OCIS codes: 160.2750, 190.4400, 190.5330 



Optical diffraction can be compensated in waveg- 
uides or when solitons form. The waves must have 
a specific size, shape, and specific intensity, for the 
solitons. Recent experiments in electromagnet ically- 
induced-transparency and nonlinear optics have identi- 
fied a situation in which diffraction is compensated, as 
for solitons, but without the scale-dependent constraints, 
a phenomenon termed diffraction cancellation [1,2]. In 
photorefractive crystals this has been demonstrated in a 
regime of so-called "scale-free optics" supported by the 
formation of a dipolar glass, a phenomenon typical of dis- 
ordered ferroelectric cystals that host glass-forming po- 
lar nanoregions (PNRs) when rapidly cooled [3,4]. Scale- 
free optics opens the way to a number of enticing effects, 
such as wavelength-insensitive propagation [5] and scale- 
free spatial instability [6]. Moreover, for its range of va- 
lidity [7], the underlying model predicts a regime that 
supports subwavelength beam propagation [8]. 

Experimentally, rapid cooling gives rise to scale-free 
optics only when the cooling rate is above a thresh- 
old value and the final equilibration temperature coin- 
cides with the Curie point Tq (in cooling), where the 
paralectric-ferroelectric phase transition occurs [2]. In 
fact, the scale-free model requires the photorefractive 
diffusion field to cancel diffraction, a phenomenon that 
entails an enhanced static dielectric response as only ob- 
served at Tq [9]. The key is that the dielectric anomaly 
is rendered accessible by the PNR-driven glass which 
suppresses long range order and the optical scattering 
typical of the equilibrium phase-transition [10,11]. 

In PNR-driven dipolar-glasses thermal history can ac- 
tually shift the crystal Tq with its associated anomalous 
enhancement of response. This property gives rise to a 
marked thermal hysteresis that delimits what is termed 
the crystal glassy phase (see Fig. la) [4,12-14]. In this 
letter we program the temperature at which scale-free 
optics occurs by exploiting specific sequences of heating 
and quenching stages. In our results, we are able to span 
the greater part of the crystal glassy phase. 

To grasp the physical underpinnings of our experi- 
ments, we recall that photorefraction leads to a diffu- 



sive nonlinearity [15,16], which profoundly alters beam 
propagation, in that diffraction is governed by an effec- 
tive refractive index n e R [2,5,8] 

n eff = n /(l - (L/Xf), (1) 

where no is the unperturbed refractive index, A the wave- 
length and L = 47rnQ£ov^XPNR(^£^Y#)- Here g is the 
effective quadratic electro-optic coefficient, Xpnr is the 
effective history-dependent low-frequency dielectric sus- 
ceptibility of the dipolar glass, Kb the Boltzmann con- 
stant, T the crystal equilibration temperature (i.e. the 
temperature measured at a given istant) and q is the 
charge of the photoexcited carriers. Eq. (1) is valid for 
L < A. As L — >> A, n e ff ^> no and diffraction is can- 
celled, the scale-free regime, independently of beam size 
and intensity. The condition L = A requires Xpnr ~ 10 5 
(which imposes T = Tq for rapid cooling schemes). 

We use a 6 x3 x2.5mm sample of ferroelectric 
Li-enriched Cu-doped potassium-tantalate niobate 
(KTN:Li). In Fig. la we report the quasi-static dielec- 
tric response e r versus T — Tc (Tc = 14.5° C) obtained 
from capacitance measurements at slow heating and 
cooling rates (\a\ = \AT/At\ ~ 0.01 °C/s) [14]. The 
marked hysteresis in e r (T — Tc) (shaded region in 
Fig. 1) signals a non-ergodic phase several degrees 
above and below Tc, a well-known complex dielectric 
response [17, 18] that is analogous to that observed in 
other relaxors [4] . The Cu doping gives it a greenish tint 
and a strong photorefraction [2]. 

Thermal preparation required to morph Xpnr is 
achieved through a computer-controlled Peltier-junction 
that fixes the temperature versus time T = T(t) sched- 
ule (the sequence of quenching and heating preparatory 
stages). Once the thermal preparation has been com- 
pleted, we shine laser light into the sample and detect 
diffraction. We focus a A = 633 nm TEMqo linearly 
polarized (in the horizontal ^-direction) z-propagating 
beam onto the input facet and analyze the input and 
output beam intensity distribution Full-Width-at-Half- 
Maximum (FWHM) using an imaging system and a 
CCD camera. The condition L = A is found when in- 
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put and output FWHM measurements coincide, i.e., the 
diffraction of micron-sized beams ceases, independent of 
size and peak intensity. 

We first carry out experiments using the standard con- 
stant a thermal preparation (rapid cooling). As reported 
previously [2], only a small region in the experimentally 
available parameter space (T—Tc,®) = (T— Tc, T) gives 
rise to L ~ A, tagged in Fig. lb as "Glassy". The region 
is evidently smaller compared to the thermal range in 
which dielectric hysteresis is detected (shaded region). 
A typical result for parameters external to the glassy re- 
gion is shown in Fig. 2b. Here, the rapid cooling process 
is characterized by (T -T c ,a) = (2.3°C, -0.12°C/s). 
As seen in the transverse intensity distribution images, 
the input intensity distribution with a FWHM of the 
Ax ~ Ay ~ 12 /im (see Fig. 2a) diffracts and spreads to 
Ax ~ Ay ~ 30 fim after the l z = 3 mm propagation in 
the sample (Fig. 2b). As expected, diffraction cancella- 
tion is not observed. Comparing diffraction to expected 
Gaussian beam diffraction (Ax out ~ 33 /im, given the 
crystal no=2.38 and A = 633 nm), we find L/X < 0.1 
(n e ff ~ no). We consider paraxial diffraction where a 
standard launch/detection scheme is sufficient (see [2]). 

We next repeat the experiment reaching the same final 
temperature (T — Tc = 2.3 °C), but through a different 
thermal preparation, characterized by a non-monotonic 
thermal path. The set of successive points occupied by 
the sample in the (T — Tc, a) plane at the different in- 
stants of time is the curve S in Fig. lb. This represen- 
tation allows to appreciate how, in terms of the rapid 
cooling experiments, the crystal is actually exposed in se- 
quence to all three "phases" , the paraelectric, the glassy, 
and the ferroelectric one. After this preparation, the 
output beam intensity distribution, reported in Fig. 2c, 
manifests no diffraction, signalling L/X ~ 1 even though 
the final equilibration temperature is external to orig- 
inal Glassy region (but still contained in the shaded 
one). We thus conclude that the spiralling trajectory S 
in the (T— Tc, a) plane is able to produce the giant Xpnr 
which is necessary to obtain L/X ~ 1 even though the 
final temperature T ^ Tc. For comparative purposes, in 
Fig. lb we also plot the R curve associated to the rapid 
cooling process (of Fig. 2b) [19]. 

We next extend the use of these "spiralling" trajecto- 
ries to obtain a scale- free response (i.e., excited glassy be- 
havior) throughout the shaded region in Fig. 1. Results 
reported in Fig. 3 demonstrate the activation of scale- 
free propagation regime at tunable temperatures beyond 
the Glassy region. We underline that our results are as- 
sociated to aging and are hence transient in time. For a 
beam with a peak intensity of I p ~ 6 W/cm 2 , the scale- 
free phenomenon occurs after an exposure of 20 s and 
begins decaying after approximately 130 s. Remarkably, 
once the aging has washed out the scale- free regime, the 
response can be rejuvenated at the same or at a different 
temperature. Repeating the full cycle, i.e., washing out 
previous photorefractive charge displacement (resetting 
of the dipolar glass), the same trajectory will lead to the 
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Fig. 1. (a) KTN:Li static dielectric constant e r for slow 
cooling (red curve) and slow heating (green curve) as a 
function of T. (b) Representation of the parameter plane 
(T — Tc,a) of the sample preparation. Both regions of 
paraelectric and ferroelectric behavior, which are unable 
to support scale-free propagation through rapid cooling, 
are indicated, together with the glassy region (dashed 
line indicates the approximate phase separation). The 
shaded region indicates the range of hysteresis, while 
the "Glassy" region indicates the temperature-cooling 
rate range where standard constant a preparation (rapid 
cooling) can give rise to to L c± A (see text). The red 
(R) and blue (S) curves represent respectively the rapid- 
cooling (a = constant) and the spiralling trajectories 
with initial (T* = 24.5 °C) and final (T f = 16.8 °C, dot- 
dashed line) temperatures that refer to the experimental 
data in Fig. 2. 



same optical response. 

In terms of soft-matter physics, our experiments 
demonstrate a method to outdo the strong temperature 
selectivity required in rapid cooling experiments through 
the complexity of the dipolar glass in the region of the 
non-ergodic phase, identified by thermal hysteresis in 
the dielectric susceptibility [3,4]. Our technique, based 
on the use of non-monotonic thermal preparation, ap- 
pears analogous to the cross-over or Kovacs effects in 
soft-matter [12, 13]. Triggered by non-monotonic ther- 
mal preparation, cross-over is known to lead to effects 
that are absent through rapid cooling [14]. 

In summary, we have demonstrated how to use cross- 
over and aging effects to tune the specific temperature 
at which scale-free optics is obtained in a disordered fer- 
roelectric. The result is a first direct example of pro- 
grammable optics, i.e. of how complexity can be used 
to achieve a desired response without resorting to mate- 
rial engineering. Loosening temperature selectivity sim- 



2 





-60 -30 30 60 

x(nm) 



0.0 



-30 30 
y (urn) 



• 

' 






60 -30 

x(nm) 


30 60 


J 





(b) 



60 -30 30 
y(nm) 




Fig. 2. Intensity distributions and x and y profiles of 
the beam in the conditions of rapid cooling (R) and 
spiralling (S) of Fig. lb. (a) Input facet of the sam- 
ple (Ax ~ Ay ~ 12 jam); (b) output facet after the 
rapid-cooling process (curve R, L/X <C 1); (c) output 
facet after the spiralling thermal trajectory (curve S, 
L = A). Note that (b) manifests a quasi-linear diffrac- 
tion (Ax ~ Ay ~ 33 /im), while the intensity profile in 
(c) is devoid of spreading (Ax c± Ay ~ 12 jam). 
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Fig. 3. Programming scale-free response at temperatures 
T ^T c . (Top right) Representation in the (T-T c ,a) = 
(T — Tc,T) parameter plane of the different conditions 
of scale-free optics reported. Activation, testified by the 
diffraction-free output intensity distribution (top) re- 
ported in the series of images Si_5, is achieved through 
spiralling thermal trajectories (bottom). Programmed 
temperatures T ~ 14.8 °C (Si), 15.7 °C (S 2 ), 16.8 °C 
(S 3 ), 17.7 °C (S 4 ), 18.7 °C (S 5 ). In the parameter plane 
representation, a is the mean value along the trajectory. 
The points lie in a temperature range of more than 4 °C. 



plifies further experimental endeavors, such as the all- 
important search for sub- micron beam propagation [7,8], 
where spatial solitons intrinsically break down [20]. 
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